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We demonstrate how, from the point of view of energy flow through an open system, rain is 
analogous to many other relaxational processes in Nature such as earthquakes. By identifying rain 
events as the basic entities of the phenomenon, we show that the number density of rain events per 
year is inversely proportional to the released water column raised to the power 1.4. This is the rain- 
equivalent of the Gutenberg-Richter law for earthquakes. The event durations and the waiting times 
between events are also characterised by scaling regions, where no typical time scale exists. The 
Hurst exponent of the rain intensity signal H = 0.76 > 0.5. It is valid in the temporal range from 
minutes up to the full duration of the signal of half a year. All of our findings are consistent with the 
concept of self-organised criticality, which refers to the tendency of slowly driven non-equilibrium 
systems towards a state of scale free behaviour. 



I. INTRODUCTION 

Water is a precondition for human survival and civili- 
sation. For this reason, measurements on water resources 
have been recorded for several centuries. A time series 
from the Roda gauge at the Nile reaches back to the year 
622 AD §. The main focus of analysis has historically 
been on statistics yielding a reliable estimate for the rain- 
fall during the growth season. The most obvious question 
to ask is in this context: How much does it rain, on aver- 
age, in the relevant months? Questions of this type can 
be answered using long time series without high temporal 
resolution, and a measurement of relatively low sensitiv- 
ity may be sufficient. Entirely different levels of resolu- 
tion and precision are needed in order to penetrate fur- 
ther into the complexity of precipitation processes. One 
might want to know just how reliable - or in fact how 
meaningful - is an estimate of future rainfall based on 
averages from the past. Of course, one would ultimately 
like to understand the processes that make a cloud release 
its water. Questions of this kind point to the statistical 
properties of rain events rather than temporal averages. 

In Sec. H we discuss the new type of radar measure- 
ment on which our analyses are based. A time-series 
of high precision rain rates with one minute resolution 
was obtained. Section [II is subsectioncd and introduces 
the various measures we apply to the time series. In 
Sec. [II A we introduce the fundamental concept of rain 



addition, we note a non-trivial relation between the du- 
ration and the size of events. In Sec. [II D we define 



events as sequences of non-zero rain rates, which enables 
comparison with many other relaxational processes en- 
dowed with an event-like structure Q. Equipped with 
this concept we investigate the statistical properties of 
event sizes. Over at least three orders of magnitude of 
event sizes the number density is consistent with a de- 
caying power law implying that there is no typical event 
size. We find that the most frequent small events are 
considerably below the typical sens itivity thr eshold of 
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and [II C 



standard rain gauges g. In Sec. 
consider the event durations and the waiting times be- 
tween successive events. The number densities of both 
event durations and waiting times follow power laws. In 



the binary signal in time of either rain or no rain and 
relate the probability distribution of waiting times to 
the fractal dimension of this signal. It is then specu- 
lated that the physical reason for the lower breakdown 
of the observed fractal regime at a time-scale of the or- 
der of 10 min may be set by the time it takes for cloud 
droplets to grow into raindrops. The upper end of the 
scaling region coincides with the time scale given by pass- 
ing frontal weather systems. In Sec. [II E we determine 
the Hurst exponent of the rain signal as 0.76, spanning 
four orders of magnitude r e [10 min, 1/2 year], extend- 
ing Hurst's result from the Nile gauge at Roda, which 
is vahd for t e [1 year, 1080 years]. Section [V estab- 
lishes a close analogy between the observed characteris- 
tics and other relaxational processes such as earthquakes 
and avalanches in granular media. Finally in Sec. ^ we 
conclude that the framework of self-organised criticality 
may serve as a useful working paradigm when dealing 
with rain. 



II. MEASUREMENT 

The recent developments in remote sensing techniques 
have opened entirely new opportunities to rain analysis. 
By using a radar rather than a common water gathering 
device, the limits on rain measurements due to evapora- 
tion, sensitivity threshold, averaging times and accessi- 
bility can be pushed considerably Q. 

The data we used refer to a height range of 50 m at 
250 m above sea level and have been collected from Jan- 
uary to July 1999 with the Micro Rain Radar MRR-2, 
developed by METEK Q . The radar is operated by the 
Max-Planck- Institute for Meteorology, Hamburg, in Ger- 
many at the Baltic coast in Zingst (54.43°N 12.67°E) un- 
der the Precipitation and Evaporation Project (PEP) in 
BALTEX . The retrieval of the rain rate is based on a 
Doppler spectrum analysis described by Atlas . At ver- 
tical incidence, the fall velocity of a droplet can be iden- 
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tified with the Doppler shift. The friction force acting 
on a faUing drop increases approximately proportional 
to its surface, but the gravitational force increases pro- 
portionally to its volume. Therefore, in the atmosphere, 
larger drops fall faster than smaller ones, and spectral 
bins can be attributed to corresponding drop sizes. For 
a given drop size, scattering cross sections can be calcu- 
lated by Mie theory ||^. Droplets are approximated by 
ellipsoids with known axis ratio ||^ . The influence of the 
changing air density with height is considered according 
to Beard llQ], and standard atmospheric conditions are 
assumed pT| . Attenuation of radar waves by droplets 
is accounted for using the observed droplet spectrum of 
the lowest range gate to estimate attenuation for the fol- 
lowing one. For higher gates all observed and corrected 
spectra of lower layers are taken into account. Thus, from 
the Doppler spectrum alone one can infer the number of 
drops rii of any desired volume Vi as well as their fall 
velocities Vi. The rain rate can be calculated instanta- 
neously as q(t) = ^^riiViVi. In the time series we in- 
vestigated, the continuous measurement is averaged over 
one- minute intervals, leading to one minute temporal res- 
olution. When the signal due to rain becomes indistin- 
guishable from the background noise at the receiver, the 
rain rate is defined as zero. Under the pertinent condi- 
tions, the calculated rain rate was typically qmin — 0.005 
mm/h, when this happened. What is measured at this 
sensitivity threshold would probably more sensibly be la- 
belled the turbulent motion of drizzle through the atmo- 
sphere, rather than rain. Instead of asking whether rain 
can be detected, a question that arises now is what we 
actually mean by rain. To achieve this level of precision, 
a conventional pluviometer would have to be able to de- 
tect a water column of 83.3 nm "rain" spread out over one 
minute. For comparison, the diameter of a single water 
molecule is about 0.3 nm. Thinking in terms of accumu- 
lated water column in such a rain gauge and given one- 
minute averaging time, one would come to the conclusion 
that the smallest detectable rain event corresponds to a 
minute during which on average every second, a film four 
molecules thick, drifts down towards the ground. This, 
of course, would be impossible to detect. The MRR-2, 
however, employs a method that is not based on water 
hitting or passing through an area of a few decimeters 
across. Given the 50 m height range starting at 250 m 
above the radar with 2° opening angle, it measures what 
happens to the liquid water in a volume of the order of 
1000 m'^, and one must bear that in mind to understand 
the minimum values calculated above. Of course, events 
at the radar's sensitivity threshold are far from being 
detectable by any water-collecting pluviometer and sim- 
ilarly far away from what we associate with the word 
"rain" . Nonetheless, we will consider any minute with 
derived q{t) > qmin as "rain", and conversely, only if the 
radar fails to detect any net downward motion of water 
through the air, we will speak of "no rain" . We will come 
back to this point in Sec. HI A. Especially for small rain 



The quantitative retrieval is restricted to rain. The re- 
flection spectra of snow and hail look very different from 
those of liquid water and can be identified. In this case 
the method fails to calculate correct water masses. The 
latest version of the instrument recognises non-rain pre- 
cipitation by an internal algorithm. The rain intensity 
data used were calculated from measurements performed 
while the development of the instrument was still ongo- 
ing, and hence the raw data had to be checked manually. 



III. DATA ANALYSIS 

The months January and February contain several in- 
stances of snow at our chosen measuring height of 250 
m. By far the largest snow disturbance was observed 
on March 6, from 3:49 am until 11:38 pm. The Doppler 
spectra reveal that 250 m altitude was inside the melting 
layer, and the water column resulting from interpreting 
the event as rain would have been 279 mm, which is of 
the order of the usual rainfall of eight weeks. In June and 
July, five very short periods of extremely high calculated 
rain rates were found (see Fig. |^). The Doppler spectra 
indicate two different types of drops with fall velocities 
at w 4 m/s and « 9 m/s. Comparison with meteorologic 
records shows that around these times, thunderstorms 
with hail or extreme rainfall may have caused the radar 
to malfunction. As in the case of snow disturbances, 
data gathered durin g thes e periods were excluded from 



HIE. The results in Sec. 



the analyse s in Sec. Ill A and 
IIIB, IIIC, and [IIEl| , however, refer to the entire data 
set since the value of the rain intensity is irrelevant here. 
To make sure that our results are not an artifact of the 
observed anomalies, all analyses were also performed on 
the certainly clean months of April and May. No differ- 
ences to the previously obtained results were observed. 
Due to the high resolution, not even the ranges of valid- 
ity were significantly affected. 
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rates the employed method is extremely powerful. 
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FIG. 1. The rain rate [mm/h], averaged over one minute, 
plotted versus the time of occurrence [min since 01-01-99, 
0:00]. The five high peaks at the right-hand side of the fig- 
ure, shown as dashed hues, are a result of malfunctioning 
during extreme weather conditions in thunderstorms. They 
correspond to 75 out of 266, 611 minutes, which is so small a 
fraction that none of our results would be significantly altered 
by including them. 

A. Event Sizes 

Previous work focused on rainfall during fixed time in- 
tervals and on the statistical properties of such fluctuat- 
ing rain intensities. Other studies addressed distributions 
of wet and dry spells (see e.g. (l2)). The fundamental 
novelty of the present study is to acknowledge the event- 
like structure of rain j^] . Events are defined as a sequence 
of non-zero rain rates, and their size M = q(t)l^t^ 
with Ai = 1 min, is the accumulated water column dur- 
ing the event. The intervals of zero rain rate between 
events are called drought periods. Our perspective is 
motivated by work on other Natural phenomena, such 
as earthquakes, where one is mainly interested in the 
events. While the entire agricultural sector depends on 
a sufficient amount of rain spread out over the months 
of the growth season, no one depends on the average 
seasonal flow of energy through the earth's crust. Due 
to this difference in anthropogenic interest, the two per- 
spectives have been used almost entirely separately in 
the respective fields. Owing to the precision and high 
temporal resolution of the data, an investigation into the 
fine structure of rain events was made possible, and the 
results §,re strikingly clear. 
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FIG. 2. The number density N(M) of rain events ver- 
sus the event size M (open circles) on a double logarithmic 
scale. Events are collected in bins of exponentially increasing 
widths. The horizontal position of a data point corresponds 
to the geometric mean of the beginning and the end of a bin. 
The vertical position is the number of events in that bin di- 
vided by the bin size. To facilitate comparison with future 
work, we rescaled the number of events to annual values by 
dividing by the fraction of a whole year during which the 
data were collected. The experimental data are consistent 
with a power law N{M) oc M~'^'^' ,tm ~ 1.4 (solid line) over 
at least three orders of magnitude, M e [Mmin,Mmax] with 
Mmin ~ 5 * 10^'^ mm and Mmax ~ 35 mm. The arrow indi- 
cates the typical sensitivity threshold of a conventional high 
precision tipping bucket rain gauge. Not only can we see that 
the radar technique is roughly 10,000 times more precise but 
also that a considerable fraction of rain events must be missed 
with conventional methods. 

Figure shows the number density of rain events per 
year N(Mj versus event size M on a double logarithmic 
plot. In a scaling regime Mmin < M < Mmax extending 
over at least three orders of magnitude, the distribution 
follows the simple power law 

N{M) cx A/-^", TM ~ 1-4. (1) 

This implies that a typical scale of events does not ex- 
ist, and scale invariance prevails. In the scaling region, if 
we compare the frequency of events of size M to that of 
events of size kM we obtain the same fraction, indepen- 
dent of M. From Eq. (|l|), it follows that: 

N{M)/N{k * M) = k^", M e [M™„, M™,,]. (2) 

But Fig. 1^ contains even more information. For events 
smaller than Mmin ~ 5 * 10"'^ mm the power law breaks 
down. This is indicative of a different physical process 
being responsible for events in this realm. Within the 
scaling regime, events of all sizes look alike when com- 
pared to others. Hence there is no reason to assume dif- 
ferent physical origins. We will later motivate the sugges- 
tion that this common origin is sudden relaxation, bursts 
of intermediately stored energy leaving the atmosphere. 
Where the power law breaks down, a different type of 
process sets in. Events smaller than Mmin might be due 
chiefly to the inner dynamics of the atmosphere. Virga, 
drizzle that evaporates before reaching the ground, is dif- 
ficult to interpret from the event perspective. Drizzle 
can form at the lower edge of clouds but immediately 
re-evaporate. Commonly the distinction between cloud 
droplets and rain drops is made in terms of diameter. 
When the droplet diameter surpasses 0.1 mm one speaks 
of rain drops. This definition reflects a physical sepa- 
ration apparent from a gap in the drop size distribution 
around 0.1 mm diameter |1^. Fringes of virga, half cloud 
and half rain, may be the explanation for events smaller 
than Mmin in Fig- Indicated with an arrow in Fig. ^is 
the typical sensitivity threshold 0.1 mm of high precision 
tipping bucket rain gauges. The value 0.1 mm is widely 
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used as the definition of zero precipitation Given 
that our interpretation of the breakdown of the power 
law is correct, and every rain event with Af > 5 * 10"'^ 
mm is actual rain, it is evident that measurements with 
today's standard precision simply don't see a consider- 
able fraction of the rain events. Questions regarding the 
fine structure of rain and the actual physical processes 
involved are then hard to address. With the radar mea- 
surement on the other hand, all rain seems to be captured 
and we can choose a suitable limit (Mmm) below which 
events are ascribed to a different physical process. 

To ascertain that we are capturing the entire physically 
relevant range of the observables of the process of rain, 
it is evidently necessary to use observational techniques 
enabling us to see beyond the physical limits of rain. Re- 
sults from investigations that do not fulfill this require- 
ment cannot be conclusive and must be treated with care- 
ful scepticism. The present study suggests a reasonable 
maximum sensitivity threshold of around 5 * 10^'^ mm, 
which is one twentieth of the commonly used threshold. 

Assuming Eq. (|l]) , we can easily calculate the number 
N{M > All) of expected events exceeding a given mass 
Ml. 

1 

N{M > Ml) oc / M-^'-'dM = Afi"^*^+\ (3) 

TAf - 1 

It follows that 

N{M > M2) = N{M >Mi){-^] . (4) 

Since we know how many events there are with M > 
Ml = Mjnin, Eq- (|) can be used to estimate N{M > 
M2), where M2 > M^in- We observed 10 events in the 
largest non-empty bin ranging from 17 mm to 35 mm, 
but from extrapolating the power law as outlined above, 
we would another 10 in the following bin ranging from 
35 mm to 70 mm. In total we would expect to see 38 
events larger than the largest event that was actually 
observed. We therefore conclude that the sudden upper 
cutoff apparent in Fig. ^ is not due to the limited time 
of observation but rather reflects a physical limit to the 
process of rain at the given location. We define M^ax as 
the largest event in the data set, a downpour of M„iax = 
35 mm of rain. 



B. Event Duration 

The number density of events versus event duration Te 
was found to approximate to a power law, see Fig. ^ 

N{Te)'^T^^^, rij«1.6. (5) 
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FIG. 3. The measured distribution of event durations 
(open circles). The data are consistent with a power law decay 
(solid line). The exponent of te ~ 1.6 is larger than the ex- 
ponent for the event size distribution, implying a non-trivial 
relation between event duration and the rain rate during the 
event. 

To highlight the implications of this result we consider 
the simplest form of a precipitation model. Naively, one 
might divide the number of minutes with measured rain 
by the total number of minutes observed, and simply 
use this fraction as the rain probability Prain in every 
minute. About 8% of the minutes we observed con- 
tain rain. Therefore, the probability for two successive 
rain minutes would be plain — 0.0064, and for 5 succes- 
sive minutes, it would already be negligible. Any model 
based on independent events produces characteristic time 
scales. In this case, one of the characteristic time scales 
would be the typical rain duration tr . The probability for 
a rain event of duration Te is given by p{Te) = P^ain- 
This can be re- written as p{Te) = e""^^/*'', where tr is 
the characteristic rain duration, which hardly any events 
will surpass. It follows that t,. = ^ in[p ) ~ 24 s. 

But the measured distribution is qualitatively differ- 
ent. Not only does the power law like number density 
allow for events longer than 1000 min, but no typical du- 
ration is found at all. We do not observe an exponential 
distribution of any kind. 

The exponent in the power law relating the duration 
of events to their frequency is different from that for the 
event sizes. This implies a non-trivial relationship be- 
tween the duration and the average rainrate during an 
event. If we could simply assume an average rain rate, 
equal for all rain events, the size would be proportional to 
the duration and the distributions would have the same 
exponent. Apparently, longer rain events are more in- 
tense. 

The statistical support for a difference between the ex- 
ponents of event size and duration is not very strong but 
the results shown in Fig. |^ reinforce this conjecture. Fig- 
ure ^ shows the average event size plotted versus event 
duration, and an exponent slightly greater than 1 is ob- 
served. This is a crude and somewhat forced measure to 
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apply but it yields results that are qualitatively consis- 
tent with Figs. H and ||, since if the average event size 
increased proportionally to the duration, the observed 
exponent would be 1. 
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FIG. 4. The event size versus event duration. The dots rep- 
resent single events of the corresponding duration. For each 
duration, the average of the single dots is evaluated (dashed 
line). Up to about 200 min the average event size {M)(t) 
increases to a good approximation like T^'^ (solid line). An 
exponent greater than 1 is consistent with Figs. ^ and ^ 
Hence, on average the rain rate was greater for long events. 



C. Drought Duration 

In Fig. H the probability distribution of drought du- 
rations N(Td) is shown to follow a power law: 



D, 



(6) 



No cut-offs were apparent. The power law is a good 
approximation from the minimum (1 minute) all the way 
to the maximum (two weeks) of observed drought du- 
rations. The only observed deviation at droughts of 
around one day length is due to the daily meteorolog- 
ical cycle. As for the event durations, this behaviour 
clearly implies correlation. We can define the drought 
probability as Pd rouah t = 1 — Prain- Hence all the ar- 
guments in sec. [II B apply for drought durations too. 
With pdrought replacing praim the typical drought dura- 
tion td = —7—7 — N ~ 12 min. The dashed line in Fig. 

a was generated by another method. Instead of treating 
minutes as the independent entities, we determine the 
rate at which rain events start by dividing the number 
of minutes by the number of rain events. This treatment 
takes into account the clustering of zero rain rates on the 
time axis, i.e. the persistence of droughts, but it cannot 
pay tribute to the dependencies which produce the real 
power-law behaviour. If arithmetic rather than exponen- 
tial behaviour persists for more than two weeks then rain 
rates at times ti and t2 two weeks apart still cannot be 
treated as independent. 
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FIG. 5. The open circles show the number density N{Td) 
of drought periods per year versus the drought duration Td- 
The solid line represents a power law approximation, with 
exponent td ~ 1.4, to the observed distribution. The arrow 
indicates one day, around which a deviation form pure power 
law behaviour can be observed. This is due to the daily mete- 
orological cycle. For comparison, a Poisson process, yielding 
an exponential distribution of waiting times, was fitted to the 
data (dashed fine). The rate of events A is defined as the 
total number of observed rain events Ntotai divided by the 
total time of observation ttotai- The number of events is then 
normalised to annual values. The Poisson process would give 
a number density 7V(10000 min) = 1.2 * 10~^* (not shown). 
Clearly, the observed values are incompatible with such an 
uncorrelated process. 

Adding the persistence of rain to that of droughts, the 
signal q{t) can be modelled with a two-state Markov pro- 
cess. One then defines transition probabilities from rain 
to drought and drought to rain, consistently with the 
fraction of total rain and drought times. In this case, 
typical drought and rain durations can be chosen. Persis- 
tence is now accounted for, but the probability for observ- 
ing drought or event durations above these characteristic 
time scales would still decay exponentially, while remain- 
ing constant for shorter events and droug hts (see Fig. |). 
This is incompatible with the observed distributions for 
both the drought durations and the event durations. The 
following section will strengthen this result further. 



D. Fractal Dimension 

A fractal (see e.g. [|4|) is a structure di splayi ng scale 
invariance of the type mentioned in Sec. 
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ing into a fractal with a factor of b and then re-scaling 
the coordinate system with a factor of h'^f , where df 
is called the fractal dimension, leaves the structure un- 
changed. Fractals often occur naturally, in which case 
the unchanged property is usually a statistical one. The 
rain data are from one fixed location but they span a 
long period of time. We define a binary signal - either 
rain or drought - and determine its fractal dimension in 
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time, using the box counting method: Different lengths I 
of time intervals (boxes) are used to cover the rainy sec- 
tions on the time axis. The number of boxes n{l) needed 
to cover the rain is proportional to . 

The results are displayed in Fig. ^. In the double- 
logarithmic plot wc find an S-shaped curve. The dashed 
lines indicate two regimes with trivial slope, dj — 1, and 
the solid line a non-trivial regime where df « 0.55. 
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FIG. 6. The number of time intervals (boxes) needed to 
cover the rain versus the box size. The fractal dimension df 
is minus the slope of this function in a double logarithmic plot, 
d/ ~ 0.55 in a scaling regime spanning 2 orders of magnitude. 
Outside the scaling regime it assumes the trivial value 1. 

Consider again the simple model with the two state 
Markov chain. As long as the box size is below the 
typical rain duration, the number of boxes needed to 
cover the rain decreases trivially; they are used to fill 
the compact space of the rain events. When the typi- 
cal rain duration is passed, each rain event is essentially 
covered by one box and the number of boxes remains 
constant. As the box size approaches the typical length 
for droughts, the entire duration of the measurement is 
filled, and the number of boxes begins to decrease in the 
trivial fashion again. The only way to obtain a non- 
trivial fractal dimension is to have - in a sense - typical 
droughts at all time scales. This amounts, of course, to 
having no typical drought duration at all. Mathemat- 
ically, this scale-freedom is represented by the power- 
law distribution of drought durations. The number of 
boxes needed to cover the rain signal will be the true 
rain duration plus the time spanned by droughts that 
are shorter than the box size (these will be overlooked), 
all divided by the box size. Hence, apart from a constant, 
representing 8% of the total time, the time Tc spanned 
by the boxes to cover the rain will increase with I as 

T, = JoN{Td) * TodTo oc /q T^^'*' * TodTo, which 
is implied by Fig. |^. Evaluating the integral we have 
Tc cx f-^^. The number of boxes needed is T^/l = r^-^^ 
In this sense a fractal relation like the one shown in Fig. 
^ could be a consequence of a power-law distribution of 
drought durations like in Fig. ||. It is the scale freedom 
that stretches the transition between the regime where 



the temporal resolution suffices to register the droughts 
and the regime where it does not. The values we measure 
suggest that there is more to the rain - no rain signal than 
only the power law of intcroccurrcnce times. Deducing 
the fractal dimension from the drought distribution only, 
we would expect a value of 0.42. But we observe 0.55, 
and the difference appears to be significant. 

The scaling regime extends from a lower limit around 
10 minutes to an upper breakdown near 3 to 4 days. 
While one might expect the fractal regime to span fur- 
ther for longer time series, the analysis of a 30 year time- 
series from Uccle suggests that the observed break- 
down is not an artifact of the shortness of our data-set. 
The authors place the cut-off at 3.5 days, which coin- 
cides with our value. Apparently, the correlation that 
gave rise to the fractal relation does not hold for longer 
than 3.5 days. Investigation of time series from Denmark 
with 1-day resolution, collected from 1876 until 2000 |l^] 
suggests that the power law for droughts does not hold 
for drought durations exceeding the upper cut-off in the 
fractal dimension. 

The explanation for the upper cutoff of the fractal 
regime may be that the typical duration of a frontal sys- 
tem moving in from the Atlantic is of the order of 3 days. 
Measured rain parameters will not belong to the same 
frontal system if the measurements are temporally sepa- 
rated by significantly more than three days. The lower 
breakdown around 10 min could not be observed in the 
Uccle time series since there the temporal resolution was 
only 10 minutes. We are still unsure as to how to in- 
terpret this lower breakdown. Clearly there must be a 
lower breakdown somewhere, and we expect it to occur 
where the particular kind of correlation that gave rise to 
the fractality on hourly to daily time scales ceases. The 
lower breakdown indicates that 10 min is a time scale 
which is special, and it must be related to a physical 
process. The microphysical processes of coagulation that 
trigger a cloud to release its water content take place on 
this time scale. Starting with typical small cloud droplets 
with radius r w 10"'^ mm, the process of stochastic col- 
lection during which small droplets merge to form rain 
drops of appreciable fall velocity takes roughly 10 — 30 
min under typical warm cloud conditions [|l6| . It is pos- 
sible that coagulation starts at a certain level inside a 
cloud and then pauses at that level before a single drop 
has left the cloud. If it then starts again, it is possible 
that on the ground we observe two layers of rain sepa- 
rated by a vertical distance corresponding to up to « 10 
min fall time. While this seems like two different events, 
from the cloud's perspective it is really only one, since 
the process of releasing water did not stop at any mo- 
ment everywhere within the cloud. Effects of motion of 
the cloud relative to the ground are not included in these 
considerations. It is unlikely that the 10 minute time 
scale is a result of the employed measurement technique. 
The radar only picks up drops with appreciable fall ve- 
locity, V > 0.5m/s. Thus the limit on the time resolution, 
given by the height of the scattering volume, is 50 m/0.5 
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(m/s) = 100 s, for the slowest drops. 



R{t)/S{t)(xt 



(10) 



E. Hurst Exponent 

In an attempt to determine the necessary size of a wa- 
ter reservoir that would never empty nor overflow, Hurst 
considered an incoming signal q(t), corresponding to 
the rain intensity in our case, that causes the level of a 
reservoir to rise or fall. Using our data, the deviation 
from the average water level in an imaginary reservoir 
would be 

t 

X(t,r)==^(q(t)-(<z).)At, (7) 
where A t = 1 min and 

(9)r = 7E'?W- (8) 
t=l 

The quantity -{q(t))r in Eq. (Q) can be thought of as 
an average outflux from the reservoir and insures that 
for any period t the water level starts and ends at zero. 
Overall trends during the interval t are thus eliminated. 
Figure shows X{t,T) as derived from the data set in 
Fig. [l[ The range of water levels the reservoir has to 
allow for is then given by 

i?(r) = max X{t,T) - min X{t,T). (9) 

l<t<T l<t<T 
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FIG. 7. Water level X{t, r) in mm in an imaginary reservoir 
for r=266,611 min, as derived from Fig. 0. During drought 
periods, a constant, slow decrease in the water level is ob- 
served, whilst during rain events the water level increases 
rapidly. The necessary size of a sufficiently large reservoir 
is given by the range indicated by a dashed line. 

Hurst determined the dimensionless ratio R{t)/S{t) 
as a function of r, where S{t) is the standard deviation 
of the influx q(t) in the period t. It can be shown that 
if q{t) is any random signal with finite variance p^ , this 
ratio increases as 



where = 1/2 is called the Hurst exponent. Hurst's 
analysis on data from the Roda gauge at the Nile, how- 
ever, yielded a different exponent of w 0.77. This 
unexpected result is commonly interpreted as a sign of 
persistence in the signal, or even as correlation. The ex- 
ponent obtained from performing the same analysis on 
our data is « 0.76 (see Fig. |). 




T [min] 

FIG. 8. The dimensionless ratio R{t)/S{t) versus r (open 
circles) shown on a double logarithmic scale. The slope of 
the fitted straight line (solid) reveals the anomalous Hurst 
exponent: R{t)/ S{t) oc with H ~ 0.76. The data deviate 
from the power law fit below r ~ 10 min in the lower limit, 
but no upper limit of the relation is observed. 

Hence, the fluctuating rain rate alone produces an 
anomalous Hurst exponent, and the result obtained by 
Hurst is valid not only for the range of 1 year < r < 
1080 years that he considered but in fact also holds for 
T = a few minutes to r = 1/2 year. Interestingly, the 
Hurst exponent deviates from this relation for t < 10 
min, which is of the same order as the observed short- 
time trivial regime of the fractal dimension. 

To understand more precisely what is actually mea- 
sured by the Hurst exponent, we applied the same 
method to a signal generated by swapping events and 
droughts at random. We kept the sizes and durations 
of rain events and droughts as determined from the real 
data and pasted them one after the other in random or- 
der. The Hurst exponent was mot altered by this proce- 
dure. In this sense it is not a measure of correlation since 
it is not affected by the order in which events occur. In 
exactly what sense it measures persistence is part of our 
ongoing research. 

IV. CONTEXT 

Self-organised criticality offers the appropriate frame 
work for dealing with relaxational process with burst-like 
behaviour whose statistics are determined by scale invari- 
ant power laws ||lj,|l^ . The term self-organised criticality 
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refers to the tendency of many systems driven by an en- 
ergy input at a slow and constant rate to enter states 
characterised by scale free behaviour. The statistics of 
the system then resemble those of a closed system near 
the critical point of a phase transition. 

A well-known example of such a process is the energy 
flow through the lithosphcre, including the outermost 
crust of Earth. Tectonic plates are driven at a slow rate 
by currents in the asthenosphere, the liquid part below 
the lithosphere, which transports heat by currents. The 
energy transferred to the plates is intermediately stored 
in the form of tension until it is finally released in an 
earthquake. Earthquake statistics follow the Gutenberg- 
Richter power law that relates the seismic moment, a 
measure of the released energy to the probability of such 
an earthquake 1 20 1 . Rather than a typical size with expo- 
nentially fewer larger than smaller quakes, scale invariant 
behaviour is observed. 

Given the right shape of grains, rice piles exhibit self- 
organised critical behaviour ||2l|. Potential energy is 
added to the system by dropping rice grains onto the 
pile at a slow and constant rate. Due to the friction be- 
tween individual grains, the pile builds up until its slope 
reaches a critical value. In this critical state, within the 
limits set by the system size, avalanches of all sizes are 
observed. During an avalanche, potential energy that was 
intermediately stored in the system is suddenly released 
in the form of heat. The distribution of energy release is 
once more a power law. 

Experiments on droplet avalanches show that self- 
organised criticality need not be restricted to granular 
media |2|]. Thresholds that enable the accumulation of 
energy before the release are given by surface tension and 
interface friction with other media. 

Rain showers share many of the features of the above- 
mentioned systems (see Tab. |). Two well separated time 
scales are present: The durations of drought periods, dur- 
ing which water evaporates, range up to months, while 
rain events take place on a much shorter time scale. The 
atmosphere receives a slow and constant energy input 
from the Sun's radiation. The absorbed energy evap- 
orates water from the surface, which is intermediately 
stored in the atmosphere. Note the analogy between liq- 
uid water in the atmosphere, tension in tectonic plates, 
and mass above the ground level in a granular pile. Dur- 
ing a rain shower, the water mass that was slowly evapo- 
rated into the atmosphere, is suddenly released, and with 
it the original evaporation energy, i.e. the condensation 
energy. The power law observed for the size distribution 
of rain events is perfectly equivalent to the Gutenberg- 
Richter law in earthquake statistics. Just as plates of the 
lithosphere don't move smoothly along, there is no con- 
stant light rainfall balancing the evaporated water mass 
immediately at every moment in time. Rain events are 
relaxations in the sky. 



V. CONCLUSION 

New insight into the working of rain can be gained by 
defining rain events, which can be regarded as energy re- 
laxations similar to earthquakes or avalanches. Taking 
this perspective, scale-free power-law behaviour is found 
to govern the statistics of rain over a wide range of time- 
and event size scales. Where clear deviations from the 
observed power laws and fractal dimensions are found, 
the limits and peculiarities of the underlying dynamical 
system become apparent, and physical insight is gained. 
Rainfall time series cannot be reproduced by conven- 
tional methods of probability theory. To enable anything 
more than an explicit reproduction of the fractal proper- 
ties, a deeper understanding of self-organising processes 
leading to fractality must be sought. Our findings suggest 
that rain is an excellent example of a self-organised criti- 
cal process. Rain is a ubiquitous phenomenon, and data 
collection is relatively easy. It is therefore well suited for 
work on self-organised criticality. For our purposes, the 
remote sensing technique employed by the MRR-2 has 
proved extremely powerful. The radar is capable of even 
higher temporal resolution than 1 min, limited only by 
the finite height of the scattering volume, and achieves 
outstanding precision in the low- intensity limit. Com- 
parison with data from other measuring sites, especially 
from warmer regions without snow and regions with more 
periodic climate would be useful in order to answer ques- 
tions regarding the universality of the observed features. 
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TABLE L Rain events are analogous to a variety of relaxational processes in nature. The two best known examples of such 
processes, earthquakes and avalanches in granular media, are summarised above. 



System 


Crust of Earth 


Granular Pile 


Atmosphere 


Energy Source 


Currents in asthenosphere 


Addition of grains 


Sun 


Energy Storage 


Tension 


Gravitational potential 


Evaporated water 


Threshold 


Friction 


Friction 


Saturation 


Release of Energy 


Earthquake 


Avalanche 


Rain event 
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